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Detailed results of heat capacity and magnetization mea- 
surements are reported on a single crystal sample of the 
spin S—l quasi-one-dimensional Heisenberg antiferromagnet 
Ni(C 5 Hi4N2)2N 3 (PF 6 ), abbreviated NDMAP. From these re- 
sults, we constructed the magnetic field (H) versus tempera- 
ture phase diagram of this compound which exhibits the quan- 
tum disordered Haldane, field-induced long-range-ordered 
(LRO) and thermally disordered paramagnetic phases. The 
phase boundary curve separating the paramagnetic and LRO 
phases is anisotropic; the increase of the Neel temperature 
(In) with H along the a and b axes is more rapid than that 
along the c axis. We calculated Tn as a function of H by 
taking into account the inter chain coupling as the form of a 
mean-field. The calculation of the staggered susceptibility of 
the S=l one-dimensional antiferromagnet with an easy-plane 
anisotropy shows quite different behavior in the different field 
directions, resulting in the anisotropic phase boundary curve. 
A good agreement is obtained between theory and experiment 
using the exchange and anisotropy constants obtained from 
the neutron scattering experiment. 

75.30.Kz, 75.fO.Jm, 75.40.Cx, 75.50.Ee 



I. INTRODUCTION 

Many lower-dimensional antifcrromagncts show a sin- 
glet ground state and an energy gap in their excitation 
spectrum originating from quantum many body effects. 
The first indication of the quantum energy gap came 
out in 1983. Haldane [EJ conjectured that the excitation 
spectrum of one-dimensional (ID) Heisenberg antiferro- 
magnets (HAFs) with integer spin quantum number (5) 
has an energy gap between the ground and first excited 
states, while the corresponding system with half-odd in- 
teger 5 has no energy gap. This conjecture has been 
tested both theoretically g and experimentally ||. It 
is now generally accepted that the quantum energy gap 
(Haldane gap) does exist in ID HAFs with integer 5. 

In real materials, there always exists an interaction be- 
tween spin chains. We call this class of materials a quasi- 
one-dimensional (Q1D) magnet. Usually Q1D antiferro- 
magnets exhibit a long-range ordering (LRO) at finite 
temperature due to the inter-chain coupling (J 7 ) The 
compound Ni(C 2 H 8 N 2 ) 2 N0 2 (C104), abbreviated NENP, 
is a typical example of 5=1 Q1D HAFs J5||| and does 



not show any indication of LRO down to 300 fiK M, 
whereas the S=l Q1D HAF CsNiCl 3 shows LRO at -4.7 
K H . The different behavior between the two compounds 
comes from the difference in the strength of J'. How ro- 
bust is the Haldane phase in a Q1D HAF against per- 
turbation? Sakai and Takahashi || studied theoretically 
the ground state properties of an 5=1 Q1D HAF with 
a single- ion anisotropy of the form -D5 2 . They showed 
that the Haldane disordered phase exists for z J'/ J<0.05 
in a rather wide range of D values, where z is the number 
of neighboring chains and J is the intra-chain exchange 
interaction. 

As was demonstrated experimentally |l(| , strong mag- 
netic fields destroy the Haldane gap and the system re- 
covers magnetism. Then, we expect a magnetic ordering 
to occur in an 5=1 Q1D HAF under high fields and at 
low temperatures. The heat capacity measurement made 
on NENP in magnetic fields did not reveal any indication 
of LRO |llj . One of the reasons for this is that an exter- 
nal magnetic field applied along the chain axis of NENP 
induces a staggered field on the Ni 2+ sites because of the 
presence of two crystallographically inequivalent sites for 
Ni 2+ [ 12) . This staggered field causes a small energy gap 
near the transition field from the disordered to the mag- 
netized states |T^.[l4|] and thus prevents the occurrence of 
LRO at low temperatures. 

We have reported the observation of a field-induced 
LRO in the 5=1 Q1D HAF Ni(C 5 Hi4N2) 2 N 3 (C104), ab- 
breviated NDMAZ ffjL a nd in Ni(C 5 Hi 4 N 2 ) 2 N 3 (PF 6 ), 
abbreviated NDMAP ]16|. Since there is only one site 
for Ni 2+ in NDMAZ and MAMAP, these compounds are 
ideal for studying the field- induced LRO. From a heat 
capacity (C p ) measurement on a single crystal sample of 
NDMAZ, we have observed an anomaly at about 0.6 K 
and 12 T EB] which indicated that a magnetic ordering 
occurred there. Because of limitations in our calorime- 
ter, we were unable to follow how the position of the 
anomaly in C p changes with temperature (T) and mag- 
netic field (H). We then tried to synthesize the 5=1 
Q1D HAF, NDMAP, with a weaker J than NDMAZ so 
that LRO is expected to be induced at a lower field. 
We have measured T and H dependence of C p in a sin- 
gle crystal sample of NDMAP and constructed the H-T 
phase diagram pq| . We found an anisotropy in the phase 
boundary curve separating the paramagnetic and LRO 
phases, and we presented a qualitative interpretation for 
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this anisotropy [[16| . Electron spin resonance (ESR) mea- 
surements made on a single crystal of NDMAP gave fur- 
ther evidence for the existence of the field-induced LRO 
phase and of the anisotropy in the phase boundary curve. 


In this paper we report detailed results of heat capac- 
ity and magnetization measurements made on NDMAP. 
We also report a theoretical analysis of the H-T phase 
diagram. 

The format of this paper is as follows. In Sec. II we 
present the relevant background information and details. 
The experimental results are given in Sec. III. In Sec. IV, 
a theoretical consideration is given on the H-T phase di- 
agram. The last section (Sec.V) is devoted to discussion 
and conclusions. 



II. PRELIMINARY DETAILS 

The compound Ni(C 5 Hi 4 N2}2N 3 (PF 6 ) (NDMAP) has 
the orthorhombic structure |18fl with the space group 
Pnmn shown in Fig. 1. The lattice parameters are, a— 
18.046 A, 6= 8.7050 A and c= 6.139 A @. The struc- 
ture consists of Ni(C5Hi4N2)2N3 chains along the c axis. 
These chains are well separated from each other by PFq 
molecules. All the Ni 2+ sites in a chain are equivalent, 
i.e., only one site exists for Ni 2+ . 

From the analysis of the magnetic susceptibility data, 
the following values are obtained p6| ; J/ks— 30.0 K, 
D/J= 0.3, .911= 2.10 and g±= 2.17, where g\\ and g± are 
the g values parallel and perpendicular to the chain c 
axis, respectively. Neutron inelastic scattering measure- 
ments were done on single crystals of deuterated NDMAP 
at T= 1.4 K JhJ. From the analysis of the data, the 
values of exchange and anisotropy parameters are deter- 
mined to be, J= 2.28 meV (= 26.5 K), J' x = 3.5xl0~ 4 
meV (= 4.1xl0~ 3 K), J' y = 1.8xl0~ 3 meV (= 2.1xl0~ 2 
K) and D= 0.70 meV (= 8.1 K), where J' x and J' y are 
the inter-chain exchange interactions along the a and b 
axes, respectively. 

The single crystals of NDMAP used in this study 
were grown from an aqueous solution of NaN 3 , 
Ni(N03)2-6H20 and l,3-diamino-2,2-dimcthylpropane. 
After filtration, KPFg was added to the solution. Well 
shaped blue single crystals up to 5mm x 5mm x 20mm 
were obtained after several weeks. Fully deuterated sin- 
gle crystals of NDMAP (NDMAP-d 2 s) were grown from 
a D2O solution of the same ingredients except that a 
deuterated l,3-diamino-2,2-dimethylpropane was used. 

The single crystals thus obtained were checked by a 
four circle X-ray diffractometer. We confirmed that the 
lattice parameters of our crystal are almost identical with 
those reported before Q . We measured C p of NDMAP 
and NDMAP-d 2 s and found that the H-T phase dia- 
grams of the two systems are essentially the same. This 
means that the magnetic parameters in NDMAP and 
NDMAP-d 28 are not different. 



Heat capacity measurements were performed with a 
MagLab HC microcalorimeter (Oxford Instruments, UK). 
The temperature and magnetic field ranges accessi- 
ble with this calorimeter are, 0.45 K<T< 200 K and 
0<if<12 T. Magnetization measurements were done 
with a MagLab VSM vibrating-sample-magnetometer 
(Oxford Instruments, UK). The temperature and mag- 
netic field ranges available with this magnetometer are, 
1.5 K<T<300 K and 0<ff<12 T. 

III. EXPERIMENTAL RESULTS 

A. Heat Capacity 

Figure 2 shows the temperature dependence of the heat 
capacity of NDMAP, including the contribution of the 
lattice measured in zero field. The data are well ex- 
pressed by the following equation, 

C p = aT + bT 3 (1) 

in the temperature range between 2 and 5 K with a= 
0.109 and 6= 0.00653. We use hereafter this bT 3 term to 
subtract the contribution of lattice heat capacity as has 
been done by many authors. 

We show in Figs. 3(a)-(c) the temperature dependence 
of magnetic heat capacity (C m ) of NDMAP, after sub- 
tracting the lattice heat capacity, in magnetic fields ap- 
plied parallel to the a, 6, and c axes, respectively. In 
all field directions, we see an anomaly in C m at finite 
fields above a critical value. This anomaly signals that a 
magnetic ordering occurs there. 

One of the advantages of our calorimeter is that a field 
dependent C p can be measured under a constant temper- 
ature |pb]| . Strictly speaking, we need to change temper- 
ature to measure C p . However, the temperature incre- 
ment necessary for the measurement is 0.5 - 1 % of the 
temperature we set so that temperature change during 
the measurement may be considered as small. Figures 
4(a)-(c) show such "field scan" data measured at several 
temperatures for the magnetic field directions parallel to 
the a, 6, and c axes, respectively. In addition to the 
sharp peak at the field denoted by -ffLRd a broad fea- 
ture is seen around the field named as H c . Here, -Hlro 
means the position of H at which the field-induced LRO 
occurs for a given T. 

Combining all the information obtained from the heat 
capacity measurements, both of the "temperature scan" 
and "field scan" procedures, we present the H-T phase 
diagram of NDMAP in Fig. 5. In addition to the 
anisotropic phase boundary (curve (A)) separating the 
disordered and LRO phases |16| , we have another bound- 
ary (curve (B)) separating the Haldane and the disor- 
dered phases which is also anisotropic. The two curves 
(A) and (B) seem to merge at a finite H when extrapo- 
lated to T= K, for respective field directions. 
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B. Magnetization 

Figures 6(a)-(c) show the temperature dependence of 
susceptibility (magnetization divided by applied mag- 
netic field, M/H) in NDMAP measured in magnetic 
fields applied parallel to the a, b, and c axes, respec- 
tively. The behavior of the susceptibility at H— 1 T is 
reminiscent of that taken at a much lower field (H= 0.01 
T) Jl(| ; a broad peak around 35 K and a steep decrease in 
susceptibility with decreasing temperature below about 
20 K. On increasing H, M/H does not extrapolate to 
zero with T^O. This behavior of M/H is similar to the 
one observed in an 5=1 ID HAF {§ and indicates that 
a transition from the gapped to a gapless phase occurs 
at a higher field. On increasing H further, M/H shows 
a minimum and an up turn at low temperatures. The 
insets of Figs. 6(a)-(c) show the low temperature part 
of the data. We see in the inset of Figs. 6(a) and (b) 
that M/H becomes almost temperature independent be- 
low a temperature whose value is field dependent. The 
temperature independent susceptibility reminds us of the 
perpendicular susceptibility (x±) of an anisotropic anti- 
ferromagnet below the Neel temperature (In). We show 
below that we are actually observing x± m this com- 
pound. Because the sign of the single-ion anisotropy term 
(DS^) is positive in this compound, spins in the ordered 
phase are expected to lie in a plane perpendicular to the 
c axis (the quantization axis of the D term is taken par- 
allel to the c axis). The anisotropy in the c plane of the 
form E(Sl — S?) is very small fPfl ]. Therefore, when H 
is applied along the a or b axes, spins point perpendic- 
ularly to H in the c plane keeping an antifcrromagnetic 
arrangement, thus giving xx- 

We plot in Fig. 5 the transition points obtained from 
the M/H data shown in Fig. 6(b). These points are 
defined as the temperatures where d(M/H)/dT shows 
a minimum for a given H. We see that the transition 
points determined from the heat capacity and magneti- 
zation measurements agree well with each other. 

IV. THEORETICAL CONSIDERATION ON THE 
PHASE DIAGRAM 

In this section, we try to reproduce the H-T phase 
diagram observed in NDMAP, exhibiting an interesting 
behavior of Tn as a function of H: At low temperatures, 
there occurs LRO only above a certain critical field and 
Tn shows an increase with increasing fields, the rate of 
which depends on the direction of H. We focus our atten- 
tion especially on the physics behind the phase diagram. 
To this end, we adopt the mean-field approximation for 
the interchain interaction which is known to work 
quite well except for fields in the vicinity of the critical 
field. In the following, we use the energy unit J=l, so 
that t, h, f and d are renormalizcd quantities of T, H, J' 
and D, respectively. According to the mean-field theory, 



the renormalized Neel temperature £n (h) as a function of 
the renormalized field h is given by the solution satisfying 
the following equation: 

l/j , z = Xrt(*N(h);h), (2) 

where Xst{t', h) is the staggered susceptibility for the one- 
dimensional magnetic system. 

Then, we calculate the staggered susceptibility of the 
magnetic chain by means of the quantum transfer ma- 
trix method combined with the finite-temperature den- 
sity matrix renormalization group. As was mentioned in 
the previous section, Ni 2+ spins in NDMAP has an easy- 
plane anisotropy, and thus the magnetic chain is well 
described by the following Hamiltonian: 

H = J[J2i S n ■ Sn+x + d{S z n ) 2 } - ]T Q^h ■ S„], (3) 

n n 

where S„ represents the spin-1 operator at the nth site. 
Neglecting the small anisotropy in the c plane, we con- 
sider the following two cases for the field h applied (i) 
along the z(c)-axis (perpendicular to the easy plane) and 
(ii) along the y-axis (in the easy plane) . 

It is noted here that the situation is quite different be- 
tween the two cases, since the former field reserves the 
axial symmetry around the c-axis while the latter breaks 
it. The nature of quantum fluctuations and hence the 
staggered susceptibility depends crucially on the symme- 
try of the system. In the following, we consider the two 
cases separately. 

(i) The field applied perpendicularly to the easy plane 

(h = (o,o,/0) 

As was mentioned, below the critical field h c the nonmag- 
netic Haldane phase called a quantum disordered phase 
is the ground state of the system, while above it the 
so called Tomonaga-Luttingcr liquid state becomes the 
ground state. Although the former has an excitation gap 
to the triplet state, the latter has a gapless excitation 
spectrum and hence is critical. This criticality of the 
Tomonaga-Luttinger liquid is characterized by the criti- 
cal exponent n defined by the divergence of the staggered 
susceptibility at zero temperature (2^] : 

Xs t(t;h)=At-^\ (4) 

where A is a constant, which scarcely depends on h in our 
calculation. Note that in the classical system r\ = 0. We 
may safely use Eq. (4) at low temperatures well below 
the temperature at which the susceptibility is maximum 
(~35 K). In Fig. 7, ry's, estimated from the numerical 
calculations for the temperature region 0.1 < t < 1.0, are 
shown by the solid circles, each of which has an error bar 
of ±0.05. The data have been shifted in fields so that the 
theoretical critical field coincides with the one observed. 
The solid curve represents the phenomenological relation 
between r\ and h: 77 = 0.3exp{— (3(h — h c )} + 0.2 with 
13 = 0.5 and h c = 0.2. 
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Now, we reproduce the phase diagram, using Eqs. (2) 
and (4) with J = 26. 5K Q. In Fig. 8 we show the the- 
oretical result by the solid curve and the experimental 
data by the solid circles. Here, we adjusted the theoreti- 
cal curve to reproduce the experimental point Tn = 0.92 
K at 11 T. From Figs. 7 and 8, we see that the in- 
crease of Tn with H is a consequence of the decrease in 
r] with H. Remembering that rj measures the degree of 
quantum fluctuations, we can say that the quantum fluc- 
tuation out of the easy plane is reduced by H so that 
the Neel state becomes more stable. In contrast to this 
quantum system, the field dependence of the constant A 
is an only source of the field dependence of Tn in the 
classical system, being very mild. We show in Fig. 7 r\ 
estimated from the phase diagram (Fig. 5) by the open 
circles, which follows also the phenomenological relation 
but [3 = 1.0. Considering the large error bars in r\ of our 
estimations, we do not think the discrepancy in f3 so seri- 
ously. Observation of the field dependence of r\ by other 
methods is desired. 

In stronger fields, spins cant in the field direction and 
hence A decreases seriously. Thus, the phase boundary 
curve closes at the upper critical field, where the mag- 
netic moment saturates. From the phase diagram, we 
estimate the interchain coupling j'z to be 1.2 x 10~ 3 , us- 
ing the value A = 2.0. This value is favorably compared 
with that obtained from the neutron inelastic scattering 
experiments ( 2j' x +2j' y =1.9 x 10~ 3 ) |Jj. 

(ii) The field applied in the easy plane (h = (0, h, 0)) 
In this case, above h c , the ground state has the Neel or- 
der and a gap opens again in the excitation spectrum 
because of the symmetry breaking field. This situation 
is quite different from the former case. The field depen- 
dence of Tn reminds us of the soliton scenario in the 
Remembering the form of the stag- 
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gered susceptibility in the classical system, we postulate 
the following in this case: 



Xs t(t; h) = (B/t 2 )exp{a(h - h c )/t}, 



(5) 



where B is a constant. This form is confirmed by 
our numerical calculations, in the temperature region 
0.1 < t < 1.0, with the coefficient a, a little less than 
unity. Although a may represent quantum effects in the 
formation energy of soliton, we assume, for simplicity, 
the classical value 1 for it. 

Now, we again reproduce the phase diagram for this 
case using Eqs. (2) and (5). The theoretical curves are 
adjusted as before using the experimental point Tn = 2.2 
K at 12 T and the critical field value H CF = 5.7 T 
for h parallel to the a-axis, and Tn = 2.7 K at 12 T 
and Hcf = 5.4 T for h parallel to the &-axis, respec- 
tively. The different values of i?cF are due to an in-plane 
anisotropy, being neglected in this paper. The phase 
boundary curves are shown by the solid curve for h par- 
allel to the a-axis and by the dotted curve for h parallel 
to the o-axis with the corresponding experimental points, 
respectively, by the solid and the open circles in Fig. 9. 



The agreement between theory and experiment is satis- 
factory. We mention that the soliton scenario still works 
in our quantum system: The symmetry breaking field 
yields the uniaxial symmetry and hence the soliton is a 
dominant source for the fluctuations in this system. Since 
the soliton formation energy, i.e. the gap, increases with 
h — h c and hence the staggered susceptibility increases 
exponentially at low temperatures, Tn shows a rapid in- 
crease with H, the rate of which is marked contrast with 
the former case. Although the soliton scenario is effec- 
tive in our case, more sophisticated study is required to 
establish further a quantum analogue of the soliton in 
classical spin chains. 

The phase boundary curve closes also at the upper 
critical field as in the former case. We estimate the in- 
terchain coupling j'z to be 1.4 x 10~ 3 , using the value 
B = 0.6. The value j' estimated again agrees with the 
value observed [fill. 



V. DISCUSSION AND CONCLUSIONS 

In Sec. Ill, we have presented detailed results of heat 
capacity and magnetization measurements on NDMAP 
from which we have constructed the H-T phase diagram 
shown in Fig. 5. We have been successful in explaining 
theoretically the phase boundary curve separating the 
paramagnetic and LRO phases in Sec. IV using the val- 
ues of intra-chain exchange interaction and anisotropy 
constants determined from the neutron inelastic scatter- 
ing measurements |0| . The inter-chain exchange interac- 
tions estimated theoretically are close to those obtained 
from the neutron experiment JlSj ] . 

We discuss the lower field boundary separating the 
Haldane and paramagnetic phases (curve (B) in Fig. 5). 
We argue below that the anomaly in C m observed along 
this curve is due to the field dependence of the first ex- 
cited triplet |l(],|24|]. We analyzed the low temperature 
part of the magnetic heat capacity data using a two-level 
system model, with a singlet ground state and the low- 
est state of the excited triplet with an energy difference 
(A_ (H)), which gives a Schottky type anomaly. We show 
in Fig. 10, A-(H X ), A-(H y ) and A_(iT 2 ) thus obtained 
for the field directions parallel to the a(x), b(y) and c(z) 
axes, respectively. The solid curves in Fig. 10 represent 
the theoretical energy level as a function of H plj ]. Here, 
we used the value determined from the neutron inelastic 
scattering experiment |l9| for the energy gap in respec- 
tive field directions at H—>0. We see in this figure that 
the agreement between theory and experiment is satis- 
factory. 

Finally, we discuss the temperature dependence of 
M/H (Figs. 6(a)-(c)). We calculated M/H as a func- 
tion of T using the quantum transfer matrix method with 
a density-matrix renormalization group technique. We 
compare theory and experiment in Fig. 11(a) for se- 
lected values of H parallel to the a axis. Here, we used 



4 



J/fcB=26.5 K obtained from the neutron inelastic scat- 
tering study fll9|| and g=2.lA determined from the ESR 
measurement J25|| . We see a good agreement between 
theory and experiment without any adjustable parame- 
ters. We have obtained the value J//cb=30.0 K from a 
fitting of the theory with the susceptibility data at high 
temperature range above about 40 K Jl6|. Because the 
lattice parameters change with temperature p8| , it is 
not surprizing if the exchange interaction constant de- 
termined at low temperatures is different from that at 
high temperatures. Figure 11(b) shows the case when H 
is applied along the c axis. Since no ESR data are avail- 
able along this direction, we assumed the value 2.05 for 
g. The agreement between theory and experiment is not 
as good as in Fig. 11(a). Further study is necessary to 
clarify this point. 

In conclusion, we have reported detailed results of heat 
capacity and magnetization measurements on a single 
crystal sample of the S=l Q1D HAF, NDMAP. From 
these results, we constructed the H-T phase diagram 
which exhibits the quantum disordered Haldane, field- 
induced LRO and thermally disordered paramagnetic 
phases. The phase boundary curve separating the para- 
magnetic and LRO phases is anisotropic; the increase of 
Tn with H along the a and b axes is more rapid than 
that along the c axis. We calculated Tn as a function of 
H by taking into account the inter chain coupling as the 
form of a mean-field. We first evaluated numerically the 
staggered susceptibility of the S=l ID antiferromagnet 
for H applied perpendicularly to the easy plane, which 
shows, at low temperatures, a typical divergence of the 
Tomonaga-Luttinger liquid with the critical exponent rj. 
Then, we got a satisfactory agreement with the exper- 
imental results using the exchange and anisotropy con- 
stants obtained from the neutron scattering experiment 
. It is interesting to note that the transition temper- 
ature Tn is governed by the critical exponent rj of the 
Tomonaga-Luttinger liquid. On the other hand, for H 
applied in the easy plane, we invoked the soliton scenario 
and got again a satisfactory agreement with experiment. 
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FIG. 1. The crystal structure of Ni(C 5 Hi4N 2 )2N 3 (PF 6 ) ab- 
breviated NDMAP. 



FIG. 2. The temperature dependence of the heat capacity 
of NDMAP in zero external magnetic field. 

FIG. 3. The temperature dependence of the magnetic heat 
capacity of NDMAP measured at the designated fields applied 
along the (a) a, (b) 6 and (c) c axes, respectively. 



FIG. 10. The magnetic field dependence of the energy dif- 
ference (A_(H)) between the singlet ground state and the 
lowest level of the excited triplet obtained from the analysis 
of Cm- The solid curves represent the theoretical result (Ref. 
24). 

FIG. 11. Comparison between theory and experiment on 
the temperature dependence of the susceptibility (M/H) at 
finite fields for (a) H \\ a and (b) H || c. The solid lines are 
the theoretical ones discussed in the text. 



FIG. 4. The magnetic field dependence of the heat capacity 
of NDMAP measured at several temperatures. The external 
magnetic field is applied parallel to the (a) a, (b) b and (c) c 
axes, respectively. 



FIG. 5. The temperature vs. magnetic field phase diagram 
of NDMAP determined from the heat capacity measurements 
(A, O) Also shown in this figure are the transition points 
obtained from the magnetization measurements (x). LRO: 
long-range-ordered phase, H: quantum disordered Haldane 
phase. P: thermally disordered paramagnetic phase. Lines 
are a guide to eyes. 



FIG. 6. The temperature dependence of the susceptibil- 
ity (M/ H) in NDMAP measured at the designated magnetic 
fields applied parallel to the (a) a, (b) b and (c) c axes, re- 
spectively. The insets show the low temperature part of the 
data. 



FIG. 7. Field dependence of the critical exponent rj char- 
acterizing the Tomonaga-Luttinger liquid. The solid circles 
denote rj's calculated for d = 0.3 and the open circles are rj's 
estimated from the experimental phase diagram. The solid 
and the dotted curves represent the phenomenological rela- 
tions reproducing the calculated and the estimated results, 
respectively. (See text.) 



FIG. 8. The H-T phase diagram for the field applied per- 
pendicularly to the easy plane (H\\c). The solid circles repre- 
sent the experimental result while the solid curve represents 
the calculated result. 



FIG. 9. The H-T phase diagram for the field applied in the 
easy plane. The solid and open circles represent the experi- 
mental results, respectively, for the field parallel to the a axis 
and for the field parallel to the b axis, while the solid and the 
dotted curves denote the corresponding theoretical results. 
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